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O , Abstract: A reliable theoretical description of double-beta- decay processes needs a pos- 
sibility to test the involved virtual transitions against experimental data. Unfortunately, 
only the lowest virtual transition can be probed by the traditional electron- capture or (3~ - 
decay experiments. In this article we propose that calculated amplitudes for many virtual 
transitions can be probed by experiments measuring rates of ordinary muon capture ( OMC) 
to the relevant intermediate states. The first results from such experiments are expected 
I to appear soon. As an example we discuss the (3(3 decays of'^^Ge and ^^^Cd and the cor- 
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responding OMC for the Se and Cd nuclei in the framework of the proton-neutron 



O ■ QRPA with realistic interactions. It is found that the OMC observables, just like the 



2i'(3(3-decay amplitudes, strongly depend on the strength of the particle-particle part of the 
fs^ \ proton- neutron interaction. 
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The two-neutrino double beta (2i//3/3) and neutrinoless double beta (Oi//3/3) decays 



■ proceed via virtual transitions through states of the intermediate double-odd nucleus. In 
^ the case of the 2^(3(3 decay only the intermediate 1+ states are involved, whereas in the 

case of the Qv(3(3 decay all the multipole states J'^ of the intermediate nucleus are active 
(see e.g. ||T|). Theoretical calculation of these virtual transitions has been a hot subject 
since two decades and a host of different approaches have been devised to accomplish the 
task The most relevant of these models, when only the ground-state-to-ground-state 
(3(3 decays are considered, are the nuclear shell model and the proton-neutron random- 
phase approximation (pnQRPA). 

The calculation of the virtual /3-decay type transitions is a formidable task and ex- 
perimental data on the corresponding electron-capture (EC) or /?~-decay transitions from 
the intermediate states to the initial or final ground state can help in fine-tuning the 
model parameters before calculation of the double-beta-decay rates. Unfortunately, the 
EC//5~-type of measurements can only probe the virtual transition through the lowest 
intermediate J"^ state. A further experimental analysis could be done by using the (p,n) 
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or (n,p) charge-exchange reactions but the extraction of the relevant information is 
not streightforward. 

In this Letter we propose to use the ordinary muon capture (OMC) as a probing tool 
for the nuclear wave functions involved in the amplitudes of the virtual transitions of the 
j3(3 decay. The OMC process is like an EC process, except that the mass of the captured 
negative muon, is about 200 times the electron rest mass The relevant OMC 
transitions here are the ones which start from the double-even 0"^ ground state and end 
on J'^ states of the intermediate double-odd nucleus. This means that in the case of the 
j3~ (3~ decays one probes the final leg of a double-beta transition and in the case of the 
double positron decays one probes the initial leg of the double-beta transition. Due to the 
large mass energy absorbed by the final nucleus in an OMC transition, one is able to probe 
intermediate states at high excitation energies, unlike in the EC transitions. A drawback 
of the OMC method, from the theoretical point of view, is that in addition to the V-A 
part the muon recoil activates also the induced parts (weak-magnetism and pseudoscalar 
terms) of the nucleonic current. This means that the theoretical expression for the partial 
OMC rates (i.e. capture rates to a particular final state) is more complicated than the 
expressions for the virtual legs of the decay, and the OMC results as such can not be 
directly interpreted as amplitudes of the virtual transitions in the decay. Instead, the 
OMC rates can be used to test in a versatile way the many-body wave functions of the 
J'^ states of the intermediate nucleus. 

In the present work we demonstrate how one can use the partial OMC rates to probe 
the structure of the J'^ states of the intermediate nucleus. This is done by making cal- 
culations of the OMC rates and the corresponding beta-decay rates in the framework of 
the proton-neutron quasiparticle random- phase approximation (pnQRPA). As test cases 
we have chosen the final nucleus ™Se of the I3~ (3~ decay of ™Ge and the double electron- 
capture (ECEC) decaying nucleus ^°^Cd. In the former case one can test by the OMC on 
™Se the states of ''^As (this is the final virtual leg in the f3~f3~ decay of ^^Ge), the most 
relevant of which are the l"*" states and the 2~ states (the ground state of ^^As is a 2~ 
state). In the latter case one can test the l"*" states of ^°^Ag by the OMC on ^°^Cd (this 
is the initial virtual leg in the ECEC decay of ^°^Cd). Also OMC rates to other multipole 
states can be measured, but the above mentioned states (1^ and 2~) have been chosen to 
serve as demonstration of the method. In addition, some beta-decay data are available for 
the lowest intermediate excitations in these nuclei so that one can probe to some extent 
the quality of the calculations for the selected multipoles. 

The experimental possibilities to use the OMC as a probe for the virtual transitions 
of the PP decays lie on prospects to use enriched targets of nuclei hke "^^Ti, ''^Se, ^^Kr, 
96 Mo, i°°Ru, ^o^Cd, i^^Sn, ^^^Xe, etc. The need for an enriched {A, N, Z) target comes 
from the fact that even a small admixture of the [A + 1, N + 1, Z) isotope would be very 
dangerous, as the probability of the neutron-emission reaction 
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fi+{A + l,N + l,Z) ^ Uf,+ {A + 1,N + 2,Z-1) ^ + {A, N + 1,Z-1)+ neutron 
is higher by an order of magnitude, and this channel would produce the same excited 
states of the final {A, N + 1, Z — 1) nucleus as the OMC in the {A, N, Z) nucleus. If the 
enrichment is less than 80 % the measurements should be done twice, for the A and A + 1 
isotopes separately, in order to separate the neutron-emission contamination. 

Also the analysis of the gamma cascades following the capture to an excited state 
in the double-odd daughter nucleus has to be performed carefully. This kind of studies 
have been performed for the Is-Od shell nuclei in T. P. Gorringe et al. by measuring 
38 gamma-ray lines and 29 (/U,z/) transitions following negative muon capture on ^^Mg, 
^^Si, ^^P and ^^S. In the heavier, double-beta-decaying nuclei the density of the final 
states of the muon capture is larger but there already exists a host of data on the energy 
levels and their gamma feeding in the intermediate double-odd nuclei. It has to be noted 
that not all of the gamma rays need to be detected. A set of the lowest ones already 
yields information on the muon-capture rates to the lowest intermediate states and can 
thus be used to probe the virtual transitions to these states. The first experimental 
proposal concerning the measurement of the partial muon-capture rates in several double- 
beta-decay camdidates, with mass numbers ranging from A = 48 to A = 150, has been 
submitted to the Paul Scherrer Institut (PSI) in Villigen, Switzerland by the Dubna- 
Louvain-Orsay-PSI-Jyvaskyla collaboration ||^. 

The 2z//3/3 decay proceeds via the 1+ states of the intermediate double-odd nucleus 
and the corresponding expression for the inverse half-life can be written as 

[4/?(o+-o+)]"' = Gg^)T |mS|' , (1) 

where Gdgt is the integral over the phase space of the leptonic variables If. The 
nuclear matrix element M^^q^ can be written as 

, ,(2.) _ En^,n(0; II E,^(j>7 II 1+) < 1+ I 1+ > (1+ II E,^(j>7 II 0+) 

where the transition operators are the usual Gamow- Teller operators, is the 2z//3/? 
Q value, Em is the energy of the mth intermediate state, Mj is the mass energy of the 
initial nucleus, and rrie is the electron-mass energy. The overlap < 1^ | 1^ > between the 
two sets of 1+ states, which are pn-QRPA solutions based on the initial and final ground 
states, helps in matching the two branches of virtual excitations [^. 

The expressions for the Oi>(3f3 amplitudes are more involved and are given e.g. in 
im, 0. In this case all the intermediate multipoles are involved in the calculation of the 
half-life and the transition operators in the virtual amplitudes contain spherical harmonics 
coming from the multipole expansion of the neutrino-exchange potential. 
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The formalism needed for the calculation of the OMC rate is developed in Ref. p!3| . 
Here we cite the main results of Ref. [13|, according to which we can write an explicit 
formula for the capture rate as 

W = AP{aZm'f^^ { 1 q\ (3) 

with A being the mass number of the initial and final nuclei, Z the charge number of 
the initial nucleus, and is the reduced muon mass. Furthermore, a denotes the fine- 
structure constant, M the average nucleon mass, the muon mass, and q the magnitude 
of the exchanged momentum between the captured muon and the nucleus. The term P 
in equation (Rl) can be written as 
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+\/3((?vQ/2M) [^J{1 + l)/{2l + d.)M[Ql + 
+ ^l/{2l-l)M[Ql- lw-]5,-„i,„) S[S-k) 



+ d^{9YQ/M){l ( \Jl + 1W{11UI, 11 + l)M[ll + 1U+] 



+ ilW{llul, 1, / - l)M[l, I - lu-]j S[^{-k) 

+ {gA/M)M[0hp]S',^{-^)6j^ + J^igA - gp){Q/2M) 



'{I + l)/(2/ + l)M[ll + lu+] + ^Jl/{21 + l)M[ll - lu-] 

xSU~^)Sr^\' , (4) 

where the W{) symbols are the usual Racah coefficients, and the definition for /, the matrix 
elements A4.[kwu{^)] and the geometrical factors Sku{i^) and k) can be found in 

Refs. [0,0. 

Finally, we note that a rather direct bridge between the (3/3 and the OMC processes 
can be established in the limit q, Z —>■ 0, since from the definition of the OMC matrix 
elements in Ref. fO, \TM it follows that 



[101] Jf^McT, (5) 
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where Mqt is the reduced matrix element 



mgt = (i+iie^o>;iio^) (6) 

of Eq. (0). Here the 0"*" state denotes either the final (for decay) or initial (for 

/3"'"/EC/3+/EC decay) 0"*" state of the /5/3-decay process. 

In the numerical computations we have used the proton model space lp-0f-2s-ld-0g- 
Ohii/2 both for the A = 76 isobars and the A = 106 isobars. On the neutron side the same 
model space was used for the A = 76 isobars, but for the A = 106 isobars the extended 
neutron space lp-0f-2s-ld-0g-2p-lf-0h was adopted. The corresponding single-particle 
energies were obtained by using the Woods-Saxon (WS) well with the parametrization of 
0. This basis was used recently in Ref. for calculation of 0z//9/9-decay rates in the A = 
76 system. For the A = 106 system, the neutron and proton WS single-particle energies 
were adjusted near the corresponding Fermi surfaces according to data on neighboring 



proton-odd and neutron-odd nuclei. This is the same basis which was used in ||10[ for a 
succesful description of the beta and double-beta data on the A = 106 isobars. 

The nuclear hamiltonian of the two isobaric regions was obtained from the Bonn one- 
boson-exchange potential [|lll with empirical renormalization by using the phenomenolog- 



ical pairing gaps, giant Gamow- Teller resonances, and spectroscopic data on nuclei close 
to the relevant isobars. For more details the reader is referred to the articles [P, |10|- As 
in very many calculations in the recent past [jl| also the results of the present pnQRPA 
calculation in the l"*" channel for the first intermediate l"*" state depend strongly on the 
strength parameter Qpp of the particle-particle part of the proton-neutron interaction. 
This strong dependence on g^p also strongly infiuences the 2z//?/3-decay rates leading to a 



strong suppression of the 2uPP matrix element, as first discussed in |]T2[. For that reason, 
to demostrate the corresponding effects on the OMC observables, we present our analysis 
as a function of g^p for the A = 76 isobars. For the A = 106 isobars we adopt the value 
(?pp = 0.8 coming from the exhaustive spectroscopy analysis of Ref. |llO| . 

Partial muon-capture rates W were calculated for transitions ^^Se(0+s ) ''^As(l'^, 2~) 
and ^°^Cd(Ogg ) — > ^°^Ag(l+). In these calculations we have used the Goldberger-Treiman 
(PCAC) value gp/gA = 7 for the ratio of the coupling strengths of the induced pseu- 
doscalar and axial- vector parts of the charged weak current. More discussion of this mat- 
ter can be found in [Q. Furthermore, the nuclear-matter renormalized value of gA = 1-0 
was adopted for the axial-vector strength. We have plotted the correponding results in 
Figs. |l|-0 as a function of gpp. Here one can clearly see that the capture rates to the if 
and 2j" states are strongly dependent on the value of the particle-particle strength ^fpp of 
the proton-neutron force. This dependence for the higher-lying states is smaller, although 
a notable exception is the OMC to the if state in ^°^Ag, as seen in Fig. |^. In Fig. |l] 
one can see that the states I2 and if cross at about gpp = 0.85, and that for gpp values 
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beyond 1.1 one obtains very small values for the OMC rate to the first l"*" state. Beyond 
this point the recoil matrix element [01 Ip] grows unphysically large and cancels against 
the contributions coming from the other nuclear matrix elements. The unphysically large 
value of the recoil matrix element stems from the excessive growth of the ground-state 
correlations, leading to the same order of magnitude for the forward- and backward-going 
amplitudes in the pnQRPA approach. This happens close to the breaking point of the 
pnQRPA, which in this case is around the value Qpp ~ 1.2. 

Partial capture rates as well as the corresponding log ft values for f3~ transitions, for 
specific values of gpp, can be found in Table ^ The value g^p = 0.9 for the 1+ channel in 
A = 76 was chosen because it would yield a sensible log/t value for the P~ decay of the 
li state (for gpp = 1.0 the log/t is far too large). The log/t value for the unique first- 
forbidden j3~ transition from the 2j" state is practically independent of the value of gpp and 
thus the unrenormalized value of gpp = 1.0 was chosen. For the A = 106 isobars the value 
gpp = 0.8 is based on spectroscopy [TO], as discussed earlier. The obtained log/t values 



can be compared with the available data, namely log/t('^^As(2j") — >-'^^Se(0+g )) = 9.7 and 
logft(^°'^Ag{lt) ^^°^Cd(0+3 )) > 4.2 ll^. Thus, a reasonable agreement with available 
data is obtained, and the predicted OMC rates can be used for estimating the possibility of 
experimental determination of the discussed partial OMC rates. In addition to the quoted 
partial OMC rates, we list in Table |^ the total muon-capture rates for each calculated J'^ 
multipole. The evaluation of the total rates has been done for the above discussed values 
of gpp. The total OMC rates, like the partial ones, are dependent on the value of gpp. 

In Fig. ^ we present histogams of all calculated partial capture rates with the corre- 
spondig excitation energy in the final nucleus. The total capture rates of Table ^ have 
been obtained by integrating these distributions over the energy. As can be seen, the 
OMC rate distribution for 1+ capture is very different for the A = 76 and the A = 106 
isobars, the distribution peaking strongly between lOMeV and 20MeV for the A = 106 
system and below 4 MeV for the A = 76 system. The 2~ distribution concentrates on low 
energies, the 2i contribution being notable. 

In conclusion, in the present work we have studied the OMC rates as tests of the 
nuclear wave functions involved in the amplitudes of the virtual transitions of the two- 
neutrino and neutrinoless double beta decay. It was found that the OMC observables 
(partial and total capture rates to 1"*" and 2~ states) depend on the strength (gpp) of the 
particle-particle part of the proton-neutron interaction. The same was observed for the 
2i//3/3-decay amplitudes already more than a decade ago. The strongest dependence is 
associated with the decays via the first l"*" state, whereas the dependence of the 2~ states 
on gpp is much weaker than for the 1"^ states. On the experimental side, the first results 
from the measurements of the partial muon-capture rates to the double-odd intermediate 
nuclei involved in the virtual transitions of the double beta decay will be available soon. 
The first experimental proposal has been submitted to the PSI, including the A = 76 and 
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A = 106 systems discussed in this article. 
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Tables 



Table 1: log ft values of /3~ decays of three lowest states of the 1+ and 2~ multipolarities, 
together with the corresponding muon capture rates W. For A = 7Q the values ^fpp = 0.9 
(J'^ = 1+) and Qpp = 1.0 (J'^ = 2~) were used. For A = 106 the value Qpp — 0.8 was 
adopted. The ratio g-p/gx = 7 was assumed in the calculations. 
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31.17 
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Table 2: Total muon capture rates with gp/g\ — 7 for the 1"*" and 2 multipoles. 



Captures 


W^tot [10=^ 1/s] 


9pp 


^«Se(0+, ) ^ 7«As(l+) 


359 


0.9 


™So(0+, ) ^ ™As(2-) 


2020 


1.0 


io6Cd(0+J ^ io6Ag(i+) 


5385 


0.8 
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Figures 
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Figure 1: Partial capture rates ^^Se(0+g ) — >• ^^As(l^) as functions of gpp witli Qp/qa = 7. 
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Figure 2: Partial capture rates ^^Se(Ol^g ) — > ^^As(2j ) as functions of gipp with g-p/ Qk 
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Figure 3: Partial capture rates ^'^^Cd(Ogg ) — > ^°^Ag(l^) as functions of g'pp -with g-p/ ( 
7. 
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Figure 4: Partial capture rates as functions of the excitation energy in the final nucleus. 



